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The weakly nonlinear regime of transverse paramagnetic dust grain oscillations in dusty (complex) 
plasma crystals is discussed. The nonlinearity, which is related to the sheath electric/magnetic 
field(s) and to the inter-grain (electrostatic/magnetic dipole) interactions, is shown to lead to the 
generation of phase harmonics and, in the case of propagating transverse dust-lattice modes, to the 
modulational instability of the carrier wave due to self-interaction. The stability profile depends 
explicitly on the form of the electric and magnetic fields in the plasma sheath. The long term 
evolution of the modulated wave packet, which is described by a nonlinear Schrodinger-type equation 
(NLSE), may lead to propagating localized envelope structures whose exact forms are presented and 
discussed. Explicit suggestions for experimental investigations are put forward. 

PACS numbers: 52.27.Lw, 52.35.Fp, 52.25.Vy 



I. INTRODUCTION 

Dusty (or complex) plasmas (DP) have been attracting 
an increasing interest among plasma physicists for more 
than a decade due to to the appearance of many novel 
phenomena. Of particular importance is the possibility 
of the existence of new plasma configurations (states), 
due to the strong intergrain coupling, including the spon- 
taneous formation of crystal-like DP structures, when 
charged microparticles (dust grains) are trapped in the 
sheath region between the electrodes, in plasma discharge 
experiments 0, Q- Such crystals, which are generally 
horizontally arranged in a position of equilibrium which is 
levitated above the negative electrode in discharge exper- 
iments, have been shown to support longitudinal (acous- 
tic) as well as transverse (optic-like) oscillation modes 


Force equilibrium in DP crystals is "traditionally" en- 
sured by the sheath electric force which exactly balances 
gravity at the levitation height. Crystalline complex 
plasma structures have been observed in recent rf dis- 
charge experiments || in which the plasma sheath was 
embedded in an external magnetic field. Gravity com- 
pensation was thus attributed to magnetic forces, thanks 
to the paramagnetic properties of the magnetized dust 
grains. Theoretical studies then followed for the inves- 
tigation of conditions for magnetic-field-assisted crystal 
equilibria involving paramagnetic charged dust grains. 
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The role of various forces acting on paramagnetic grains 
has been discussed in Ref. [j], where magnetic forces 
(due to magnetic dipole interactions) have been shown 
to prevail over the (weaker) electric polarization forces 
(especially under experimental conditions considered in 
Ref. 0], where magnetic fields as strong as a few thou- 
sand Gauss were used). According to these considera- 
tions, the properties of transverse magnetized dust lat- 
tice (TMDL) oscillations were investigated where 
both isolated grain linear oscillatory modes and prop- 
agating linear waves were discussed. The dependence of 
the dynamics of (linear) grain oscillations on the specific 
characteristics of the (inhomogeneous) magnetic field was 
studied in detail in Ref. and the possible occurrence 
of a linear instability depending on the field profile was 
pointed out. 

It is now established that the investigation of the linear 
regime of a dynamical system only unveils part of its dy- 
namical profile. Nonlinearity may be present in DP crys- 
tal dynamics, either due to intergrain interactions or due 
to the sheath environment. In a generic manner, as the 
oscillation amplitude becomes large, nonlinearity is first 
manifested via the generation of wave harmonics leading 
to the amplitude modulation of the waves. Increasing 
displacements even further, one may come up with the 
formation of localized structures (solitons), which propa- 
gate and interact with each other in a remarkably stable 
manner, thanks to a compensation between nonlinearity 
and dispersion. The present study is devoted to an inves- 
tigation of weakly nonlinear effects, with respect to mag- 
netized transverse dust-lattice oscillations. Nonlinearity 
is explicitly shown to be related to the magnetic and/or 
electric field profile, as well as to the intergrain electro- 
static and/or magnetic (dipole-dipole) interactions. By 
using a perturbative approach, the generation of phase 
harmonics is elucidated and exact expressions are ob- 
tained for the (weak) vertical displacement of the param- 
agnetic dust particles. Once these oscillations propagate 
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in the dusty crystal as a transverse wave, the amplitude 
is shown to be potentially unstable to external pertur- 
bations, under conditions which depend explicitly on the 
magnetic field characteristics. Finally, the possible for- 
mation of localized envelope excitations is discussed, and 
the dependence of these coherent structure characteris- 
tics on the plasma parameters is pointed out. 



II. NONLINEAR SINGLE GRAIN 
OSCILLATIONS 

Let us consider the vertical motion of a charged dust 
grain (mass M and charge q), subject to an external 
static electric and magnetic field, E and B respectively, 
both in the vertical (~ z) direction. The vertical dis- 
placement 5z — z — zq from the equilibrium position Zq 
obeys the equation of motion 



F 
M 



(1) 



where we have set zq = 0. The three terms in the right- 
hand side (rhs) of Eq. (1) account for the (sum of the) 
electric and magnetic forces, viz. F = F e + F m , in the 
z— direction, the force of gravity, and the usual (Ep- 
stein) damping term, which involves the phenomenolog- 
ical damping rate v due to dust-neutral collisions. 

We shall assume a smooth, continuous variation of the 
(generally inhomogeneous) field intensities E and B, as 
well as the grain charge q (which may vary due to charg- 
ing processes |y|) near the equilibrium position zq = 0. 
Thus, we may develop 



E{z) » E 



E' z 



:E'' 



B(z) w B + B' Q z + -Bq z 2 + 



and 



q(z) sa q + q' z + 



1 // 2 

2 %Z 



where the prime denotes differentiation with respect to 
z and the subscript '0' denotes evaluation at z — zq, viz. 
Eq = E(z = z ), E' q = dE(z)/dz\ z=Zg and so forth. Ac- 
cordingly, the electric force F e — q(z)E(z) and the mag- 
netic force F m — —d(mB)/dz = —2aBdB/dz (where 
the grain magnetic moment fi is related to the grain ra- 
dius a and permeability \x via m = (/i — l)a 3 B/(/j, + 2) = 
aB 0), which are now expressed as 

F e (z) S3 q E + {q E + q' Q E ) Z 



-(q EZ + 2q' E' Q + q' ( ;E )z 2 



and 



F m (z) sd -2aB B' - 2a(B' Q 2 + B Q B" ) z 



-a(B B%' + 3B' B" )z 2 
may be combined into 



F P + F n 



<9$ 

dz 



where we have introduced the phenomenological poten- 
tial $(z) 



$(z) sa $(z ) + 
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dz 
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Ci) = {^{z)/dz\ =zo = 

{qE )^~ 1] +a(B 2 ) { Q 3) (for j = 1,2,...; the superscript 
(j) denotes order in partial differentiation) are obvious 



The definitions of $ 



(1) 

(2) 



-(qE) + a(B 2 )' Q = -q E + 2aB B' 
-(qE )' + a(B 2 )Z 



= -(q' E + q E' ) + 2a{B'l + B Q B'£) 
*(3) - -(qEoYt + a(B 2 )Z 

= -(q'jEo + 2q' E' + qoE%) + 2a(3f?X + BqB'q) , 

(3) 

and so forth. Note the force balance equation 

Mg = q E - 2aB B' , 

which is satisfied at equilibrium. 

Given the above definitions, the equation of motion |T)l 
takes the form 



z + v z + lu 2 z + Ki z 2 + K 2 z 6 



0. 



(4) 



where 



c^ = $ (2) /M, Kx = $ (3) /(2M) , K 2 = * (4) /(6M) , 

(5) 

and higher order terms are omitted. One immediately 
notices the intrinsically nonlinear character of the trans- 
verse dust oscillations due to the electric/magnetic field 
inhomogeneity. In the linear limit ($(j) = for j > 3), 
the results of Ref. are exactly recovered. 

Once the set of parameter values are determined, 
Eq. (0J can be solved, e.g. via a Lindstedt-Poincare 
method. Assuming small displacements of the form 
u(t) = Y^n=i e ™ u n( T ), where the reduced time vari- 
able r incorporates a frequency shift due to nonlinearity: 
t = Lut = (1 + eAi + e 2 A 2 + ...)uj g t = (1 + \)u) g t. One thus 
obtains the solution 



u(t) sa eA cosr+ e 2 



-±,4 2 cos2t-^L4 2 
3w 2 u> 2 



(6) 
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along with a nonlinear frequency shift 5u> — Xuj g ~ 
|^4| 2 w g . We see that the nonlinearity of the force act- 
ing on the paramagnetic dust grains typically results in 
the generation of frequency harmonics once the oscilla- 
tion amplitude becomes slightly significant in magnitude. 
The occurrence of this effect depends on the relative mag- 
nitude of the coefficients K\ and K2 in Eq. J3J|. This 
phenomenon has already been studied for an unmagne- 
tized dust crystal (viz. for B = Q) :8|; based on Eq. 
(9) therein, one deduces that ^(3) / < J ) (2) = — 1 mm -1 
and <I>(4)/$(2) = +0.42 mm~ 2 , implying K 1 /ui 2 , = -0.5 
mm -1 and K2/uJg — +0.07 mm -2 . Similar data on the 
magnetized DP crystal case, which are not yet available, 
may be deduced from appropriate experiments. This 
type of analysis may be pursued further once such feed- 
back from experiments is available. 



III. NONLINEAR TRANSVERSE 
MAGNETIZED DUST LATTICE OSCILLATIONS 

Let us now consider a layer of identical charged dust 
grains (of lattice constant ro). The Hamiltonian of such 
a chain is of the form 



^ 2 \ dt 



U{r r , 



where r„ is the position vector of the n— th grain. 
U nrn {r nm ) is a binary interaction potential function, re- 
lated to the electrostatic intergrain interaction potential 
4>{x) (typically of the Debye type, though ion flow in 
the sheath may be included for a more complete de- 
scription 0, as well as the magnetic moment of the 
ri— th and to— th grains, which are located at a distance 
fnm = I r" Tl — r m |. Even though the analytical form of 
U(r) need not be specified here, for generality, we may 
explicitly refer to the model of Refs. 0, |B| in the case 
of magnetized dust crystals [see Eq.(15) in Ref. or 
Eq. (8) in Ref. The external potential & ext (r) ac- 

counts for the external force fields in the sheath region 
(i.e. essentially $ as defined in © above); nevertheless, 
in a more sophisticated description, <& ex t may include 
the parabolic horizontal confinement potential imposed 
in experiments for stability |Toj . or the initial laser ex- 
citation triggering the oscillations in experiments (both 
neglected here). 

Considering the motion of the n— th dust grain in the 
transverse (vertical, off-plane, ~ z) direction, we have 
the equation of motion including dissipation caused by 
dust-neutral collisions 



dt 2 



8z n 



(7) 

where F ext = —d$ ext (z)/dz accounts for all external 
forces in the z— direction. 



A. Equation of motion 

Assuming small displacements from equilibrium, one 
may Taylor expand the interaction potential U(r) around 
the equilibrium intergrain distance Iro — \n — m\ro (be- 
tween Z — th order neighbors, I — 1,2,...), i.e. around 
z n « 0. Retaining only nearest-neighbor interactions 
(1 = 1), we then obtain from Eq. (6) 



J 0,T 



(2z r , 



n+1 



%n) (^ro %n— l) 



(8) 



The transverse oscillation characteristic frequency wo.t 
and the coefficient K3 are defined by 



LU ,T = {-U , (r )/(Mr )} 1 / 2 , 
h',= --^[U'(r )-r U"(r )}. 



2Mrl 



(9) 



The gap frequency Lo g and the coefficients K\ and Ki 
have been defined above. Of course, a negative/positive 
value of U'(ro)/U"(ro) is a condition ensuring stability 
of transverse dust-lattice (TDL) oscillations, i.e. a real- 
valued frequency LOo^/^g ( - /t here means • or f, respec- 
tively); cf. © above. 

Notice that the equation of motion presented in Ref. 
[f| [see Eq. (18) therein] is exactly recovered in the linear 
case (i.e. for K\ = K2 = K3 = 0), given the form of the 
interaction potential in that model. In particular, the 
frequency u ,t then becomes [BJ 



U 0.T 



Mrl 



- D ) CXP Ad 



9^ 



where Ad denotes the Debye length; one immediately 
distinguishes the contribution from the Debye potential 
(first term) from the second term, which is due to the 
magnetic moment too = aBo (see the definitions above). 
Notice, in passing, that the vertical motion equation of 
the recent nonlinear model by Ivlev et al. [Tl| is exactly 
recovered in the appropriate limit |T^ . 

In general (i.e. regardless of the particular analyti- 
cal aspects of each model) , one immediately notices that 
nonlinearity may either arise from the sheath environ- 
ment (electric/magnetic fields) or from the interactions 
between the paramagnetic dust grains. 

Adopting the standard continuum approximation, of- 
ten employed in solid state physics [l3| . we may assume 
that only small displacement variations occur between 
neighboring sites, viz. 



du 1 J 2 



1 ^dh 



1 A d A 



z n±1 = u ±r - + -r —± -r Q — + -r — ... 

, where the vertical displacement z n (t) is now expressed 
by a continuous function u = u(x,t). One may now pro- 
ceed by inserting this ansatz in the discrete equation of 
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motion (|SJ, and carefully evaluate each term. The cal- 
culation leads to a continuum equation of motion of the 
form 



Cn u 



^xx 
2 



+ c 



2 r 



U, 



- LO g U — K\U 



°12 

K 2 u 3 + K 3 rt (ul) a 



(10) 



in the following order e n+1 and so forth, we obtain the 
to— th harmonic amplitudes Un ^ at each order, along 
with a compatibility condition up to any given order. 
The equation obtained in order ~ e 1 is 



(d 2 



r 2 V 2 



Lq u\ = . 



(11) 



where terms of higher nonlinear ity have been omitted. 
We have defined the characteristic TDLW velocity co = 
^o.t^o; the subscript x denotes the partial differentia- 
tion, viz. (u%) x = 3(u x ) 2 u xx . 

The continuum equation of motion l|lUfl is a modi- 
fied, damped Boussinesq-like nonlinear equation. In the 
following, we omit the phenomenological damping term 
(which may be added a posteriori, once an explicit solu- 
tion for the displacement is obtained). 



We may assume that u\ = exp[i(kx — uit)] + c.c. = 
exp(i#) + c.c, where u>, k = 2tt/X and A denote the 
carrier wave frequency, wavenumber and wavelength, re- 
spectively; c.c. stands for the complex conjugate. The 
dispersion relation obtained from is of the form 



r 2 

fc 4 ^ r- 
12 1 



J T,0 ' 



„2 2 



1 



„2 2 



12 



(12) 



IV. AMPLITUDE MODULATION - A 
NONLINEAR SCHRODINGER EQUATION 



According to the standard reductive perturbation 
method [lj], we shall consider a small displacement of 
the form: u = eu\ + e 2 U2 + ... , where e <C 1 is a 
small parameter and solutions u(x,t) at each order are 
assumed to be a sum of m— th order harmonics, viz. 
u n = Sm=o u ""^ exp[i(fcir — ut)] (the reality condition 



,(») 



,(»r 



is understood). Time and space scales are 



accordingly expanded as 

d/dt d/dTo + ed/dTi + e 2 d/dT 2 + 



and 



d + ed!+e 2 d 2 



d/dx -> d/dXo + ed/dX 1 + e 2 d/dX 2 + ... 



V + eVi 



e 2 V, 



2ed d 1 +e 2 {df + 2d d 2 ) 



V 2 Q + 2eV Vi +e 2 (V 2 + 2V V 2 ) 



implying that 
d 2 /dt 2 -> 

9 2 /fe~ 
and 

9 4 /9a; 4 -» + 4e VgVi + 2e 2 Vg (3V 2 + 2V V 2 ) + ... . 

This reductive perturbation technique is a standard pro- 
cedure, often used in the study of the nonlinear wave 
propagation (e.g. in hydrodynamics, in nonlinear optics, 
etc.) [lit. 

We now proceed by substituting all the above series in 
HI U| l and isolating terms arising in the equation of motion 
at each order in e n . By solving for u n , then substituting 



predicting an inverse optical-mode-like behaviour (in 
agreement with previous results; in fact, this is exactly 
the continuum analogue of the discrete dispersion rela- 
tion in Ref. 0). We assume that the reality condition 
lu > is, in principle, satisfied in the region of validity of 
the continuum hypothesis k -C (ir/ro) = k cr (neverthe- 
less, note that the possibility of magnetic field-related 
instabilities, basically related to the relative magnitude 
of lo 2 and lq 2 -, , was put forward in Ref. ) . 

Let us evaluate the action of the linear operator Lq, 
defined above, on higher harmonics of the phase 6 

L e me = [(-inoj) 2 + uj 2 g + c 2 (ink) 2 + c 2 (r 2 / 12) (ink) 4 } 
xe md 

= [n 2 (n 2 -l)c 2 (r 2 /12)fc 4 -(n 2 -lV 2 ]e^ 



(13) 



where we made use of the dispersion relation (|12|l ; thus 



o = uj g , Di — and D 2 



-3a, 2 + c 2 r 2 k\ 



D, 

To order ~ e 2 , the condition of suppression of secular 
terms takes the form 



(i) 



(i) 



dTi 



dXx 



(14) 



i.e. — Ui\(,), where C — e i x ~ v gt), implying that 
the slowly-varying wave envelope travels at the negative 
group velocity v g = u>'(k) = —(1 — fc 2 rg/6) c\ k/w; this 
backward wave has been observed experimentally |16| . 



We notice that v a becomes zero at k 



a/6? 



be- 



yond which it becomes positive; nevertheless, note for 
rigor that this comment rather becomes obsolete once 
one properly takes into account the influence of the lat- 
tice discreteness on the form of the dispersion relation 
lu = uj(k); cf. Ref. [13 • 

This procedure finally leads to a solution of the form 



u(x, t) 



Ae i{kx-ujt) +cc 



5 



e 2 K, 



+0(e 3 ) . 



2\A\< 



-Do 



2i (kx—cut) , 
3 v 1 + C.C. 



(15) 



where to obeys the dispersion relation (|12|1 . 

The amplitude ^4(X, T) obeys a Nonlinear Schrddinger 
Equation (NLSE) of the form 



dA d 2 A 
i l P 

dT dX 2 



Q\A\ 2 A = 0, 



(16) 



where the 'slow' variables {X, T} are {X\ — v g Ti,T2}, 
respectively. The dispersion coefficient P, which is re- 
lated to the curvature of the phonon dispersion relation 
CE2J as P = (l/2)(d 2 cu/dk 2 ), reads 



P = - 



2cj 3 



l -2 k ' 



and the nonlinearity coefficient 



1 












2uj 


L V w 9 





3 K 2 - 3K 3 k 4 ^ 



(17) 



(18) 



is related to both the sheath electric/magnetic field and 
the intergrain coupling nonlinearity discussed above. As 
expected, (|18ll recovers exactly the previously derived ex- 
pressions (in similar studies) in the appropriate limits, i.e. 
precisely Eq. (9) in ^8| in the dispersionless limit (for 
very low k, i.e. cancelling terms in fc 4 in H12I) ). and Eq. 
(10) in [l7| (where a quasi-continuum approximation was 
adopted, i.e. a continuum envelope but discrete carrier 
wave desciption). Notice however that inter-grain inter- 
action nonlinearity was neglected in those studies (i.e. 
K3 = therein). 

In view of the analysis that follows, it may be appro- 
priate to cast Eq. I|l(jfl m a non-dimensional form, for 
physical clarity. This is readily done by introducing a 
set of appropriately chosen scales, e.g. the lattice con- 
stant ro and the inverse eigenfrequency oj~I , i.e. by scal- 
ing the (slow) variables X and T as X — > X 1 — X/ro 
and T — > T' = cj g T; the vertical displacement ampli- 
tude becomes A — > A' = A/r . The form of Eq. itPo|) 
is then exactly recovered upon substituting with: Eq. 
P -> P' = P/(io g r 2 ) and Q -> Q' = Qr 2 /tu g . The primes 
will be dropped in the following. 

a. Modulational instability. In a generic manner, a 
modulated wave whose amplitude obeys the NLS equa- 
tion l|16fl is stable/unstable to perturbations if the prod- 
uct PQ is negative/positive. To see this, one may 
first check that the NLSE is satisfied by the monochro- 
matic solution (Stokes' wave) A(X,T) = A e iQ|Ao|2T + 
c.c. The standard (linear) stability analysis then shows 
that a linear modulation with the frequency and the 
wavenumber k obeys the dispersion relation 



Q 2 (k) = P 



2 k 2 



2||A,| 2 



(19) 



which exhibits a purely growing mode if k > k ct = 
(Q/P) 1 / 2 \Aq\. The growth rate attains a maximum 
value jmax — Q |^4o| 2 - This mechanism is known as the 
Benjamin-Feir instability [l5j. For PQ < 0, the wave is 
modulationally stable, as evident from (|19fl . 

Notice that the coefficient P is negative, given the 
quasi-inverse-parabolic form of the dispersion curve cu(k) 
for low k. One therefore only needs to deduce the sign of 
Q, given by ifTHl) . in order to determine the stability pro- 
file of the TDL oscillations. In fact, given the definitions 
(13 EI) of the parameters u> g , K%, K2 and K3, one may eas- 
ily derive an expression for Q in terms of the (derivatives 
of the) potentials $(z) (for the sheath) and U{z) (for 
inter-grain interactions). The exact form of the poten- 
tial $(z) (hence the coefficients j — 1, 2, ...) may be 
obtained from ah initio calculations or by experimental 
data fitting. As mentioned above, some evidence for the 
numerical values of the coefficients $(j), yet only in the 
absence of the magnetic field, can be found in Ref. 0, 
where the dust grain potential energy was reconstructed 
from experimental data (see Eq. (9) therein); those val- 
ues (cf. comment above) seem to suggest that the value 
of Q for low wavenumbers k is positive, as may be read- 
ily checked from l)18[l. Therefore, under the experimen- 
tal conditions described in Ref. [8j , the TDL wave would 
propagate as a stable wave, for large wavelength values A. 
However, for shorter wavelengths, either the coefficient 
P = oj"(k)/2 or Q may change sign, and the TDL wave 
may thus be potentially unstable. These results should a 
priori be checked by appropriately designed experiments, 
in particular with regard to magnetically levitated DP 
crystals. 

b. Localized excitations. A final comment concerns 
the possibility of the existence of localized excitations 
related to TMDL waves. It is known that the NLSE 
(fTfi|> admits localized solutions (envelope solitons) of the 
bright (PQ > 0) or dark/grey (PQ < 0) type. These 
expressions are found by inserting the trial function A = 
Aq exp(iO) in Eq. H16(l and then separating real and 
imaginary parts in order to determine the (real) functions 
Aq(X, T) and Q(X, T). Details on the derivation of their 
analytic form can be found e.g. in Refs. 0,|2(|, so only 
the final expressions will be given in the following. Let 
us retain that this ansatz amounts to a total (vertical) 
grain displacement u(x,t) essentially equal to: 



u(x, t) = ep cos(kx — tot + 0) . 



(20) 



where p = 2Aq and the nonlinear phase-shift is O ~ e, 
since linear in {X, T} — {ex, et} (see below). 

The bright type (pulse) envelope solutions (continuum 
breathers, see Figs. HEJ, obtained for PQ > 0, are given 
by 



A = 



2P 
QL 2 



1/2 



sech 



X-v e T 
L 



exp |ii[ Ue X+(0-|)T]|, (21) 



6 



where v e is the envelope velocity; L and f2 represent the 
pulse's spatial width and oscillation frequency, respec- 
tively. In our problem, the bright-type localized envelope 
solutions may occur and propagate in the lattice if a suf- 
ficiently short wavelength is chosen, so that the product 
PQ is positive. We note that the pulse width L and the 
amplitude A® satisfy LAq — (2P/Q) 1 / 2 — const.. Let 
us point out that, when the pulse's spatial width L is 
comparable in order of magnitude to the carrier wave- 
length A, this (now highly localized) type of solution is 
similar in structure to the (odd-parity) discrete breather 
(DB) modes (or intrinsic localized modes, ILMs) recently 
widely studied in molecular chains pi) . 

For PQ < 0, we have the dark envelope soliton (hole) 
[13 (see Fig. E| 



A = ±^ tanh 



(X~v e T 



(2PQA 



!1 

2 



T] k?2) 



which represents a localized region of negative wave den- 
sity (void). The pulse width V = (2\P/Q\) 1 / 2 /A is 
inversely proportional to the amplitude A . 

For PQ < 0, we also have the grey envelope soliton 
HH (see Fig. EJ) 

A = A Q {l-d 2 sech 2 {[X -v e T]/L"}} 1/2 exp(iG), 

. (23) 

where 6 = 0(X, T) is a nonlinear phase correction to 
be determined by substituting into the NLSE (|16[) : the 
calculation yields the complex expression: 



9 



1 

2P 



v x 



-S sin 



1, 



-X-2PQA T + 9 



d tanh 



1 - d 2 sech 2 ( ^f- 



1/2 



(24) 



(see (64) in Ref. [l!|a). This localized excitation 
also represents a localized region of negative wave den- 
sity; 9o is a constant phase; S denotes the product 
S = sign(F) x sign(i> e — Vb). Again, the pulse width 
L' = (IP/QI) 1 / 2 /(dAo) is inversely proportional to the 
amplitude Aq, and now also depends on an independent 
real parameter d, which regulates the modulation depth; 
d is given by: d 2 = 1 + (v 2 - V$)/(2PQA ) < 1. V is 
an independent real constant which satisfies (see details 
in Ref. pjj): V -^/2\PQ\A < v e < V + y/2\PQ\ A . 
This excitation represents a localized region of negative 
wave density (a void), with finite amplitude at X = 0. 
For d — 1 (thus Vq — v e ), one recovers the dark enve- 
lope soliton presented above, which is characterized by a 
vanishing amplitude at X = 0. 

We should admit, for rigor, that the latter excitations 
(of dark/grey type), yet apparently privileged in the con- 
tinuum limit (where PQ < for low fc), are rather 
physically irrelevant in our (infinite chain) model, since 



they correspond to an infinite energy stored in the lat- 
tice. Nevertheless, their existence locally in a finite-sized 
chain may be considered (and possibly confirmed) cither 
numerically or experimentally. 

It may be stressed that the grain displacement cor- 
responding to the (envelope) excitations presented here 
is intrinsically different in form (and obeys differ- 
ent physics) from the pulse-like small-amplitude local- 
ized structures (solitons) typically found via Korteweg- 
DeVries (KdV) theories; see Refs. mm for a critical 
comparison. 



V. CONCLUSIONS 

The present study has been devoted to an investigation 
of amplitude modulation effects associated with either 
isolated transverse dust-grain oscillations or propagating 
transverse dust-lattice waves in DP crystals of param- 
agnetic charged dust particles, embedded in an external 
magnetic field. We have shown that nonlinearity comes 
into play once the oscillation amplitude slightly departs 
from the weak-displacement (linear) regime. This non- 
linearity, which is related to the electric and/magnetic 
field(s) in the sheath region, affects the dynamics of 
transverse dust lattice oscillations via the generation of 
phase harmonics and the potential instability of the car- 
rier wave, due to self-interaction. The latter may presum- 
ably be responsible for energy localization in the DP crys- 
tal via the formation of localized envelope excitations. 
Analytic expressions for these excitations are presented 
and briefly discussed. It should be stressed that the in- 
stability suggested here, related to self-modulation of the 
carrier wave and triggered once the amplitude becomes 
slightly important, is completely independent from the 
one pointed out in Ref. 5] , which is related to the values 
of the intrinsic parameters (the gap frequency, in partic- 
ular) involved in the physics of the problem. 
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Figure Captions 

Figure 1. 

Bright type (pulse) soliton solution {PQ > 0) of the 
NLS equation obtained for an indicative set of 

numerical values for the parameters in Eqs. i|20l 121 [I : 
e = 0.1, P = Q = 1, v e = 0.2 (i.e. ,0 = 2,4.0 = V2), 
k = 2tt/X = 2tt and L = 2: (a) the waveform, as results 
from (|2U|) ; (b) the corresponding (localized) wave energy 
(~ p 2 ) - normalized value. 



Figure 2. 

Same label and data as in Fig. except L = 0.2. 



Figure 3. 

Dark type soliton solutions of the NLS equation for 
PQ < (holes); here Q = —P = 1 and the remaining 
values are just as in Fig. ^ (a) the waveform, as results 
from (|20|) ; (b) the corresponding (localized) wave energy 
(~ p 2 ), normalized over it asymptotic value. 



Figure 4. 

Grey type soliton solutions of the NLS equation for 
PQ < 0: values identical to those in Fig. [3J in addition 
to Vq = 0.5, d = 0.9. Notice that the amplitude never 
reaches zero. 
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